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THE STRAIN OF A GRAVITATING SPHERE OF VARIABLE 
DENSITY AND ELASTICITY* 

BY 

L. M. HOSKINS 

I. Introductory. Statement of the problem 

The problem of determining the strain of a gravitating elastic sphere from 
an initial condition of hydrostatic equilibrium, under the action of small 
disturbing forces having a spherical harmonic potential, has been solved for 
the case in which the density and elastic moduli are uniform throughout the 
body. It has also been solved for cases in which the density varies with the 
distance from the center, but with the restriction that the modulus of compres- 
sion is infinite while the modulus of rigidity has a uniform value. 

This problem is of especial interest because of its relation to estimates of 
the rigidity of the earth. Such estimates are based upon the comparison of 
the actual yielding of the earth to tidal and centrifugal forces (as inferred from 
certain refined observations) with the computed yielding of an elastic sphere 
having the same size and mass as the earth. For the purpose of such a com- 
parison it is obviously desirable that the ideal sphere assumed in the computa- 
tions should agree with the actual earth-sphere as closely as possible in all the 
elements essential to the problem. Now it is quite certain that the density 
of the earth varies greatly with distance from the center, that the material of 
the earth is far from incompressible, and that the earth as a whole is very 
much more rigid than the surface rocks. Our knowledge of the actual proper- 
ties of the earth is much better represented by the assumption that all three 
of the quantities, density, modulus of compression, and modulus of rigidity 
are functions of distance from the center than by any of the more restricted 
assumptions which are essential to the above-mentioned solutions of the 
sphere problem. It is therefore desirable to solve, in as general a form as 
possible, the following problem: 

A sphere of isotropic material, in which the density and elastic moduli are 
functions of distance from the center, and which would be in hydrostatic 
equilibrium under self-gravitation alone, is strained by small disturbing forces 
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having a potential which is a spherical harmonic function of the coordinates 
and a simple harmonic function of the time; it is required to determine the 
strain. 

The analysis given in the present paper leads to a solution of this problem 
which, though not fully general, is a considerable advance in generality over 
the restricted solutions above mentioned. The solution of the general problem 
is shown to depend upon that of a pair of simultaneous ordinary differential 
equations; and it is found that these may be solved in terms of convergent 
series when certain restrictions are imposed upon the functions expressing the 
density and elastic moduli. Moreover, these restrictions are not such as to 
detract seriously from the geodynamical interest of the solution. 

The scope and order of the paper may be summarized as follows : 

I. Introductory. Statement of the problem. 

II. Establishment of fundamental equations in general form. The density 
and elastic moduli are here treated as unrestricted functions of distance from 
center, and the disturbing potential as a simple harmonic function of the time. 

III. Solution of the equations when the density is variable while the elastic 
moduli are constants.* The density function is assumed to be rational and 
integral. 

IV. Solution of equations when the density and both elastic moduli are 
variables. The assumption as to density is the same as in III, while the 
elastic moduli are represented by rational integral functions of a restricted 
form. 

V. Particular solutions. Particular forms of the functions expressing den- 
sity and elastic moduli are substituted in the formulas expressing the general 
solution given in IV, reducing these to the working formulas used in numerical 
computations. 

VI. Application to the earth. Numerical results are given for the strain 
of a sphere having the size and mass of the earth, showing the effect of different 
assumptions regarding density, elasticity, and compressibility. Besides the 
results for the static problem, a series is also given for the case in which the 
disturbing potential is a simple harmonic function of the time. 

II. Establishment of the fundamental equations 

1. Physical theory. The analysis which follows is a generalization of that 
used in a former paper dealing with a less general form of the sphere problem,! 

* The assumptions of variable density and variable elasticity are found to be independent 
of each other as regards their effect on the fundamental equations, and it seems advantageous 
to present separately the two parts of the analysis leading to the most general form of the 
equations which express the solution. 

t The strain of a gravitating, compressible elastic sphere, these Transactions, vol. 11 
(1910), pp. 203-248. 
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but with an amendment to the physical theory. This amendment leads to 
certain changes in the fundamental equations; but though physically impor- 
tant, these changes are algebraically of minor effect, leaving the same general 
method of solution still applicable.* 

The problem of the strain of an elastic solid from a condition of great 
initial stress cannot be treated by the equations of elastic deformation in the 
ordinary form, since these equations assume the proportionality of actual 
stress to strain. Some of the earlier papers on the strain of a gravitating 
sphere treated the problem as if the body were strained from an initial 
" natural " configuration in which gravitation did not act. Applied to a 
body of planetary dimensions such a solution leads to deformations far beyond 
the elastic limits of known materials; but even apart from this fact, the basal 
assumption can have no validity as applied to actual bodies. It is only on 
the assumption of incompressibility that this method of solution leads to 
valid results. If the material is incompressible the computed effect of gravita- 
tion alone is a condition of hydrostatic stress throughout the sphere, with 
zero strain; so that the determination by this method of the strain due to 
small disturbing forces combined with gravitation amounts to the same thing 
as the determination of the strain due to the disturbing forces alone from an 
initial condition of hydrostatic stress, f For a compressible material, however, 
the reasoning fails. 

It is reasonable to suppose that a solid body is elastic for small deformations 
even if initially in a condition of great stress. It seems reasonable, also, to 
apply the ordinary laws of elasticity to such a case, with this change: the 
relation " stress is proportional to strain " is replaced by the relation " in- 
crement of stress is proportional to strain."! 

* The defect in the physical theory as formerly stated was pointed out by the writer in a 
footnote to a later paper, these Transactions, vol. 11 (1910), p. 504. The same 
point was discussed by A. E. H. Love in his work Some problems of geodynamics (1911), in 
which the amended theory was applied in an able analysis of the sphere problem in the case 
of uniform density. 

t This method was employed by G. Herglotz in a paper dealing with the strain of an in- 
compressible sphere of variable density. Vber die Elastizitat der Erde bei Beriicksicktigung ihrer 
wriablen Dichte, Zeitschrift fur Mathematik und Physik, vol. 52 (1905), 
p. 275. In this paper the solution for the general case of density varying with distance from 
the center is reduced to the solution of a linear ordinary differential equation of the sixth order 
whose coefficients depend upon the law of density, and the complete solution is given for 
certain cases of that law. » When the assumption of incompressibility is introduced the analysis 
given in the present paper leads to a differential equation identical with that obtained by 
Herglotz. 

J In the application of this principle in the writer's former paper cited above, "increment 
of stress" was computed by considering a definite volume-element before and during the strain. 
This was later seen to be an error, because the displacement of an element of matter causes a 
change in the gravitational force acting on that element which is of the same order of im- 
portance as the disturbing force. The reasonable method is rather to compare the stress-condi- 
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The differential equations of strain, formed in accordance with this principle, 
involve the initial stresses; these are in general indeterminate, and the solu- 
tion of the equations is therefore indeterminate unless something is assumed 
regarding the initial stress-condition. The assumption which has commonly 
been made, and which will be made in the analysis herein given, is that the 
initial condition is that of hydrostatic equilibrium. This assumption is not 
only especially simple but is especially interesting in the geodynamical problem, 
since in the case of a body of planetary dimensions like the earth the pre- 
dominating stress throughout the greater part of the body is without doubt 
a normal pressure.* 

2. Differential equations of strain of an isotropic elastic solid from a condi- 
tion of initial stress. As a preliminary to the discussion of the problem of 
the sphere it is useful to establish the fundamental equations in as general a 
form as possible. We therefore consider first the general case of an isotropic 
elastic solid undergoing a small strain from any state of initial stress; the 
fundamental assumption being that increment of stress is connected with 
strain by the same relations which, in the ordinary theory, are assumed to 
connect actual stress with strain. 

In establishing the differential equations it is convenient to use rectangular 
coordinates; polar coordinates are however employed in the application to 
the sphere, f The notation for displacements, strains, and stresses will for 
the most part be that adopted in A. E. H. Love's treatise The Mathematical 
Theory of Elasticity, second edition; axial components of displacement are, 
.however, denoted by u x , u u , u z instead of u , v , w . 

The equations of motion for an element of the body are 

d 2 u x _ , dX x , dX v . dX z 
P -&2T = PX + ~ST + 



dft dx dy dz 

m d 2 u v ^dY x . dY v , dY t 

(1) PTZ = pY +■ 



dx dy dz ' 

d 2 u z dZ x BZ V dZ z 

P ~W = pZ + te + -dy- + ~dz- 

These equations are independent of elastic theory. 

tion of an individual element of matter before and during strain. This is the amendment to 
the physical theory referred to above. The defective theory had also been employed in a 
paper by A. E. H. Love, The gramtational stability of the earth, Philosophical Trans- 
actions of the Royal Society, London, A, vol. 207 (1907), p. 171. 

* This method of treating the geodynamical problem seems to have been first suggested 
by Lord Rayleigh. See paper On the dilataiional stability of the earth, Proceedings of 
the Royal Society, London, A, vol. 77 (1906). 

t Attention is called to the fact that after the change to polar coordinates the letters x , 
y , z are given new meanings. 
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The ordinary stress-strain relations for an isotropic elastic solid are expressed 
by the six equations 

Equations (2) apply to the case in which the body is strained from its " natu- 
ral " state, i.e., from a condition of zero stress. 

Let it now be supposed that the body is initially in equilibrium under the 
action of any bodily and surface forces, so that initial stresses exist throughout 
the material, and that u x , Uy , u t denote displacement from this initial con- 
figuration; it will be assumed that equations like (2) hold with the substitution 
for each stress-component of the increment of that stress-component due to the 
strain. 

Let SX X , dY v , • • • denote these increments of the stress-components for an 
element of material initially at (x, y , z); then the assumption is expressed 
by the equations 

SX,-\A + 2„^, 8F„-M + 2„^«, JZ. - XA + 2„^-', 
These equations give, by differentiation and addition, 



l sx * + i sx » + i sx > = ^ + ^ d i 



(4) 

+ — A 4. 2 — — 4- 
dx dx dx 



dy \ dx + dy ) + dz\dz + dx )' 



with two similar equations which may be written by cyclic permutation of 
x, y,z. 

Let the components of actual stress at the point (x, y, z) before strain be 
X X ,X V , • ••; and let X' x , X' s , • •• denote the same quantities- at the same -point 
during strain. Then the increment of stress for an individual element of 
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material initially at (x, y , z) is given by equations like the following: 

/en sv v v i dXx , dX x dX x 

(5) dX x = X x - X x + u x -^ + u y -^- + u a ,-^-; 

and these values substituted in (4) reduce the first member to the expression 

(dF, ax; dx[\_(dx 1 dXy dx,\ 

\ dx + dy + dz ) \ dx + dy + dz J 

d ( dX x dX x dX x \ d ( dX v dX y dX y \ 

(6) + rx\ U *^ + U ^ +U -d7) + Vy\ U *-te +U >'-dy- +U *-dz-) 

,a ( dx,^ dx, ax,\ 

+ dz{ U *-dx- + U "-dy- +U °-dz-)- 

But the first of equations (1) gives 

dX x dX y dX z 
~dx~ + ~di + ~dz~~ ~ pX > 

(7) , , , 

dX x dX' y dX z , ,d 2 u x 

~dx' + ~dy~ + ~dz~ =: ~ p X +P ~d?-> 

in which accented letters denote values during strain and unaccented letters 
corresponding values before strain at the same point (x, y, z) . Also 



(8) 



(dp dp dp \ 

U *dx + U »dy: + U °d-zr 



and since p, X may replace p' , X' in terms containing u x , u v , or u z , the 
expression (6) becomes 

d ( dX x dX x dX x \ d ( dX y dX y dX y \ 

(9) + dx\ Ul -dx~ +U "-di +U 'lz~ ) + dy\ U *-dx- +V *-d y - +U *-dz- ) 

, d ( dX. dX z dX.\ 

+ dz\ U '-dx- + U "~dJ + U '~dz-)' 

Equation (4), with (9) substituted for the first member, is one of three dif- 
ferential equations applicable to the general case of strain of an isotropic 
elastic solid from a condition of initial stress. The six stress-components 
whose derivatives appear in these equations are connected by the relations 
obtained from (1) (with the acceleration terms omitted) ; but the solution will 
be indeterminate unless something more is known concerning the initial 
stresses. 
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3. Simplification of equations when initial stress is hydrostatic. It will 
hereafter be assumed that the initial stress-condition is that of hydrostatic 
equilibrium, so that all terms involving tangential stress-components vanish, 
while X x = Y y = Z t . Using (1), the expression (9) may now be written 

(10) 

— ^- {pu x X + pu y Y + pu, Z) . 

4. Form assumed by equations when the body is a gravitating sphere. 
Let it now be assumed that the body is a sphere, and that p, X, and ju vary 
only with r, the distance from the center. If V denotes the gravitation 
potential before strain (a function of r only), 

mi Y- xd l Y- ydV 7- ZdV 

(11} X -~rdr> Y ~~rTr> Z rTr ' 

Hence if u T denotes radial displacement, 

dV 

(12) U x X + U y Y + U Z Z = Ur-fo) 

d_ 



(pu x X + pUyY + pu t Z) = ^ ( f> u ^ ) 

d ( dV\ 
= P d-x\ U 'Tr) 

(dp dp dp\ xdpdV 



(13) 
Also 



xdpdV 
+ Ur rdr~d7- 



dV 



Assuming now that the bodily forces during strain are derived from a potential 
V , the expression (10) reduces to the following: 

R 2w * , dV * 9 ( rr , „ , dV\~] 

Substituting this for the first member of (4), and writing E x for the terms in 
the second member which contain derivatives of X and ju> we obtain the follow- 
ing as one of the three differential equations for the strain of a gravitating 
sphere: 
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(16) 



dA 
(X + it) j^ + mV 2 u x + E x 



=>m^-i(r'-^m- 



The other two equations are obtained by cyclic permutation of x , y , z . 

The first members of these equations are the same as would be obtained if 
there were no initial stress. The equations are in fact identical! with those 
which would apply to strain from zero-stress under the action of forces whose 
axial components (per unit mass) are 



(17) 



i{y 



V + Ur 



dV 

dr 






dV\_ 
dr ) 



dV 
dy 



A, 



dV 

dz 



In other words the strain may be treated as if due to three sets of bodily forces: 

(a) Forces having potential V — V; these include the disturbing force 
and the change in the gravitational force due to the changed configuration 
of the attracting mass. 

(b) Forces having potential v* (dV/dr). 

(c) The increment of gravitational force per unit volume due to the change 
— pA in the density of the attracted element. 

5. Equations in polar coordinates for gravitating sphere. We now pass to 
polar coordinates, replacing x, y, z by r, 8, <j> as independent variables, and 
u x , Uy , u z by Ur , u 9 , u^ as component displacements, according to the following 
scheme of direction cosines. 





r 


6 


<t> 


X 


sin 9 cos i> 


cos 6 cos $ 


— sin <t> 


y 


sin 6 sin <j> 


cos 6 sin <t> 


cos <t> 


z 


cos 


— sin 6 






Each of the required equations is obtained by multiplying the three equa- 
tions like (16) by the proper direction cosines and adding. The second 
members of the resulting equations may be written by inspection (noticing 
that, to the first order of small quantities, the acceleration components are 
d i u r /dt 2 , d 2 u e /dt 2 , d 2 ujdt 2 ). The first members, aside from the terms 
obtained from E x , E v , E z , are given in treatises dealing with the case in 
which X and ju are constants.* Let the terms involving derivatives of X and \t 

* See A. E. H. Love's Treatise on the Mathematical Theory of Elasticity, second edition 
(1906), p. 138. For the values of vr r , vrg , vr^, , A , see p. 56. 



(18) 
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be E r ,E e ,E^; then the three equations in polar coordinates are 

[d*U r , dV d (-., -., dV\-\ 

X + 2 M dA 2 M f3 5w r l 

7^nT^"7L^ (rBr ' )_ ^J + * 

Here ur r , w>, ot$ are rotational components of strain. 

It remains to determine E r , E t , E+ by carrying out the process above out- 
lined. If X and n are unrestricted functions of the coordinates, the resulting 
values are 

v . d\ dju dttr , 1 d/x / 1 du r du t u t \ 
T ~ *dr + dr dr + r 30 \r dd + dr r ) 



+ 



dfl ( 1 du. 



r sin d<f> \ r sin 6 dtj> 

3/t /l dttr du$ _ 
r dd 



dr r )' 



Ad\ dn(l du r du t _u e \ 2 dfi/tir 1 du t \ 
*'- ~ *a+dr\r d0 + dr r) + rdO\r + rde) 



1 duf 1 d«» 1 dw» u* A 
(19) +rsin»30Vsinea^ + r d0 r e0tV )' 

♦ ~ r sind d<t> dr\r sin d<l> dr r ) 

1 ^M / 1 du e 1 dUj U4, \ 

+ r dO\r sind d<l> + r d$ r COt ° J 

+ 2 *f+ *. 1 «g\ 

r sin 9 3<^ \ r r r sin o<£ / 

Since in the subsequent applications of equations (18) it is to be assumed 

that X and n are functions of r only, the expressions (19) are much simplified.* 

* Considerable algebraic work is involved in the deduction of (19), but it is merely the 
routine work of transformation of coordinates and is therefore omitted. Although the applica- 
tion in the present paper is limited to the case in which X and p are independent of 6 and 4> > 
it seems desirable to record the general expressions. 
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6. The boundary conditions. The condition which will be assumed at the 
boundary surface is that the increment of stress on this surface is zero.* 
Making the same assumptions as in equations (3), but using polar coordinates, 
the increments of normal and tangential stress on a surface initially perpen- 
dicular to the radius vector are 

bu, 

m ^ (du u e 1 dUr\ 

(20) ^="Ur "7 + rWj' 



. / duj, _u* 1 du r \ 

r * ~ M V rd r + r sin 6 d<j> ) ' 



and these are to vanish when r = a. 

7. Case in which the disturbing potential is a spherical harmonic function 
of the coordinates and a simple harmonic function of the time. We now as- 
sume that the disturbing potential is proportional to Si r*, where Si is a spher- 
ical surface harmonic 'of order i and a simple harmonic function of the time. 
In this case equations (18) may be reduced to simultaneous ordinary differ- 
ential equations by the following method.t 

Assume 

(21) u r = uSi, n 9 =v^, u,=^- e ^-, 

in which u, v are functions of r only. Substituting in the known formulas 
for A, vr T , vr t , vr^ , making use of the partial differential equation 



(22) 



1 d 2 St 1 d ( . n dSi 



+ s keU sined £>) + iii + 1)Si = ' 



sin 2 d<j> 2 
and introducing the notation 

(23) »-J5— ST"— (• + !). p z = - r \Tir- U )> 
we find 

,^.x n „ „ n 2 dSi _■ dSi 

(24) A=y Si , 2*r r = 0, 2^=-^^, 2 -*-«m ; 

so that the first members of (18), exclusive of the terms E r , E e , E+, take 

* The amendment to the physical theory, referred to above, affects the boundary conditions 
as well as the differential equations. The increment of stress is zero, not for the fixed surface 
r = o , but for the material surface initially at r = a . 

t This method was used in the writer's previous papers dealing with less general cases of 
the sphere problem. Some of the details given in those papers are here omitted. See these 
Transactions, vol. 11 (1910), p. 238 et seq.; also p. 495 et seq. 
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the forms 

|_ x r r ar J sin d</> 

while from (19), 

p ^ dx „ , <*M<frA dfifu dv v\dSi 

(26) 

_ dp fu dv v\ 1 dSi 

* ~5r \r + 5r~ r/ sine a^' 



Again, since <S,- is a simple harmonic function of t , we may put 

(27) -^«-p»S„ 

in which p is a constant. Hence the acceleration components in (18) become 

< 28 ) -a?---^«*. -a? = -P 2c W ^ = -sinW- 

It remains to express the values of V and V for substitution in the second 
members of (18). 
Let the potential of the disturbing forces be 

(29) W = ^k riS <> 

c being a constant which measures the intensity of the disturbing force, and 
g the surface value of the gravitational attraction per unit mass. 
If p m is the mean density of the sphere and a the radius, 

< 3 °) V = T~n \l {' Pr*dr + fpnfr] . 

p m a |_ r J J r J 

The potential V is made up of W and the gravitation potential of the 
strained body; hence if U denotes the increment of the gravitation potential 
caused by the strain, 
(31) V'-V = W+U. 

The value of U may be expressed as the sum of two parts, one of which is 
the potential of a distribution of matter of density 

a d P 

~P*-Ur Tr 
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throughout the sphere, the other that of a surface layer of thickness equal to 
the surface value of u r . Since, by virtue of (21) and (24), 

(32) _ / a_ w g_-( /If + ll *) 5 , f 

the former of the two parts of U is 

while the latter is 
(34) 



iwySi 



i piMia (0' 



2i + 

in which the suffix G) refers to surface values. Combining (33) and (34) 
and integrating by parts the terms containing dp/dr , we find 

Wi r i r ( «. d 



U (2t + l) ft 

(35) 



A^£ p { yri+2 -d\^y 



+ r* 



in which 4iry has been replaced by 3g/p m a . The value of V — V is obtained 
by combining (29) and (35), in accordance with (31). 
From (21), (24), and (30), 

dV 3<7«S t - r dV SgySi C T , 

dr p m ar 2 J r dr p m ar 2 J 

The final result of the substitution of (21) in (18) is obtained by combining 
the foregoing partial results. The first members of the resulting equations 
are made up of (25) and (26), while the second members are obtained from 
(28), (29), (35), and (36). After cancellation of the factors Si, dSi/d$, and 
(l/sin0) (dSi/d<f>) , the three equations are independent of 0, <f>, and t, 
the second and third equations being in fact identical. The solution of the 
system (18) is thus reduced to that of two simultaneous ordinary differential 
equations with r as independent variable and u , v as the functions to be deter- 
mined. These two equations, after multiplication by r, may be written as 
follows: 

<*+*>£+«<«+i>».+'£»+*£J 



(37, _«(«_2Lj£>*)_„ 



• 2 ru, 
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< 38 > " (S?IjT.[?iI'(»' , "-i<- , ">) 4 

The substitution of (21) in (20) reduces the boundary conditions to the 
two equations 



(39) 



X,, + 2^=0, M( Jr --+-j = 0, 



to be satisfied when r = a . 

It will be shown that the solution of (37) subject to the conditions (39) 
may be obtained in the form of infinite convergent series in a comprehensive 
class of cases in which p, X, and p. are variable. The solution will first be 
given on the assumption that X and p. are constants while p is expressible as a 
rational integral function of r . It will then be shown that the effect of assum- 
ing X and tx also to be rational integral functions of r is merely to add terms to 
the main equations employed in the solution, and that when these functions 
are restricted in a certain way the series expressing u and v are still convergent. 

III. Solution of equations when the density function is rational 

AND INTEGRAL AND THE ELASTIC MODULI ARE CONSTANTS 

1. Solution when both elastic moduli are finite. Assuming X and p, to be 
constants, let 

(40) p = po(l + £ ia ; + A; 2 a; 2 + ••■), 

in which p , ki , k 2 , • • • are constants, and 

(41) x=-. 

a 

Hereafter x will replace r as independent variable, and u, v will be replaced 
by e, a, defined by the equations* 

(42) e = l «-? 

r r 

Also, since p 2 a is of the same dimensions as acceleration, let 

(43) p 2 a = ng . 



* It will later appear that, in the important case i = 2 , e and a have simple physical mean- 
i ngs. It is seen that x , e , and a (as well as y and z already defined) are abstract numbers. 
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Equations (37) and (38) now take the following forms: 

(X + 2p.)xy- + i(i + l)p.z = — \ x l I pa? dx — np m ex 2 ) , 

ax p m \ ax x ,/„ J 



(44) 



(X + 2fi)y + p. , = (if - np m ax 2 ), 

ax p m 



R = 



(45) 






2i + 

+ x* 

while the boundary conditions (39) become 

(46) Xy + 2 M ^^ = 0, M (*^ + e)=0, (whena^l). 
If we now assume 

00 00 

y = 2-t ^m % , z = 2Lj & m x , 

00 00 

e= T, A m x m , a = 23 5„af», 

m=<— 2 m=i— 2 

it is found that the coefficients C m , D m , A m , B m may be so determined as to 
satisfy identically both (44) and (46). The solution may proceed as follows. 
Notice first that A m and B m may be expressed in terms of C m and D m by 
means of (23), which may now be written 

1 d{ex?) \ d{ax?) 

(48) y = x 2 ^x~- l{% + l)a ' Z = x-lx-- e - 

These give 

C m - (m + Z)A m — i(i + l)B m , 
(49) 

D m = -A m + (m + 2)B m , 

which must hold for all values of m appearing in the assumed series (47). 
Since the series for y and z contain no terms of degree* i — 1 or % — 2 , equa- 
tions (49) require that 

(50) Ai-2 -iBi-2 = 0, 

(51) Ai.. x = £,_! = 0; 

while for all other values of m , 

* It will later appear that in the case of fluid equilibrium ( n = and n = ) , terms of 
degree i — 1 must be introduced unless ki = . This exception has no effect on the g, . :a 
method of solution. 
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A (m + 2)C m +i(i + l)D m 

m (m + 2)(m + 3) -»(t + l)' 
(52) 

r C m + (m + 3)D m 

m (m + 2)(ra + 3) -t(* + l)* 

Introducing now the notation 

(53) k' = po/p m , b = gp m a/ p. , 

we proceed to substitute the assumed values of e, a, y , z in (44) and (45). 
The expression for R takes the form 

00 

(54) R = PoT.LnX™, 

in which 

ii 

+ 3 



\ . 2t - 2 _c_ 

y i ~ 2t + 1 i_2 2&" 



(55) 

+ + 

T m -i + 4 

while for m = z+l,z" + 2, ••• 



,m+2 



3 

+ &1 ^»-3 + 2& 2 ^m-4 + 3& 3 Am-5 + ■ • • ) 
+ 3(|^ m _ 2 + ifa Amr-z + \h An-4 + •••). 

Equations (44) become 

Z[Q + 2jmC„, + ;(t+ 1 )/),,, J *- 

= bk' 2 (l + k 1 x + k i x* + ••■)lT,mL m * n ~ 3£G,(1* 
(57) + **!***+ ••O.-pE^.ar"], 

Z [ Q + 2) C m + (m + 1) Z> m ] x m 

= bk' 2 (i + hx + k 2 x>+ •)[z^^-pIX* m+2 ]; 

the summations extending to all values of m occurring in the assumed expan- 
sions (47). These equations are to be satisfied identically. 

The right-hand members contain no term of degree lower than i. The 
same is true of the left-hand members if (50) and (51) are satisfied. Hence to 
satisfy (57) identically it is necessary and sufficient to satisfy the following 
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equations f or m = i, i + 1 , i + 2 , • • • : 

Q + 2 J mC m + t(t + 1)7). = 6** [ml* + fe x (m - 1 )i m _ 1 

+ h(m-2)L m - 2 + ••■ - (c„_ 2 + fc; c m _ 3 + &; c^ + • • • ) 

(58) - p (^4m-2 + &i A m - i + & 2 ^^-4 + ■ ■ ■ ) , 

Q + 2 J C. + (m + l)7) m = W^lz* + &i L m _! + & 2 L^ 2 + • • • 



41 

~ p ( 7?„,_ 2 + &i 5 OT _3 + k 2 B m -i + 

in which k\ , k 2 , • • • are easily-determined functions of ki , k 2 , 
For m = i the two equations are identical : 



']• 



(59) Q + 2^ d + (i + l)D ( = bk*(li - ^7^.-2) ; 

while f or m = i + 1 , i + 2 , • • • , they may be solved for C m and 7) m , giving 

(m - t) (m + t + 1 ) D m = M' 2 f h 7^ + 2& 2 7^ 2 + 3& 3 7,„ 



^m— 3 



71 

+ • • • + C m _ 2 + &i C„_3 + k 2 Cm-4 + • • • + p ( ^4 m _ 2 — mB m - 2 

(60) + ^! (^ OT _ 3 - mB^-3) + & 2 ( A.-4 - mB^-i) + • • • ) , 

Q + 2^C m = - (m + 1) 7) m + 6&' 2 [z m + feZw + k 2 L m -, 



ft 

+ • • • - p (B m - 2 + ftx B m _ 3 + h 5 m _4 + • • • ) 



]■ 



These equations serve to determine any pair of coefficients CU, D m (except 
Ci, Di) in terms of coefficients of lower order. Remembering the relations 
(50), (51), (52), (55), and (56), it is seen that every C m and D m may, by succes- 
sive applications of (60), be ultimately expressed in terms of the three coeffi- 
cients Ci , D t ,Ai-2. When this has been accomplished the only non-vanishing 
terms in (57) will be those of degree i; and these will vanish when (59) is 
satisfied. Equation (59) in fact becomes one of three linear equations for 
determining C,-, 7),-, and Ai- 2 ; the other two being obtained from the surface 
conditions (46). These three equations may be written as follows: 
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(£ + 2^ C t + (i + l)Di = ^[a^n £ (°» -(rn-i + 2)A m ) 

x ( 1 + -+h + ^_ +...) 

(61) + 1^| A*+ - ± - ± A*+ ] , 

X Z C m + 2 M £ (m + 1 ) A m = 0, JU I(m5 n + .4)=0. 

m=i m=i— 2 m=i— 2 

The foregoing solution is valid for any constant values of X and y. which 
are compatible with the assumption that the strain is "small," i.e., that 
only the first powers of the displacements are of appreciable magnitude. 
The statement of the solution however requires modification in the limiting 
cases of incompressibility ( X = °o ) and fluidity ( n = ) . 

2. Case of incompressibility. If X = » the cubical expansion A must be 
, hence in applying the general solution to this case we must put y = while 
\y remains finite. Putting 

(62) -v = y', 



M 



equations (44) and (45) become 
(63) 



dy' 
x dx~ + i(i + 1)z 



gpa ( dR \ 

= -M m \ X te- nPmeX )' 



y' + J S 1=9 ^ R - n ^ <x ^' 



p- 3 

R -~2i + 



'-i[^i Xp d\ {exM)dx + xi £ p i{^) dx \ 

(64) 

+ — I px 2 dx — „cp m x\ 

X J„ I 

Without restricting the law of density, a linear differential equation for 
determining e may be obtained by eliminating y' between the two equations 
(63), substituting for z and a their values in terms of e as given by (48) with 
y = 0, and substituting for R its value given by (64). By so operating as to 
eliminate the integrals containing e, there results a linear differential equation 
of the sixth order whose coefficients depend upon p and its derivatives.* This 
equation will not be given; it may be noted, however, that the elimination of 
the integrals containing e is accomplished by operating upon (64) as indicated 

* In the case n — this equation becomes identical with that obtained by Herglotz in the 
paper cited above. His method of procedure is, however, quite different. 
Tram. Am. Math. 80c. a 
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in the following equation: 

,(PR , n dR .,.,,.„ 
(65) ** * 



" 3 (*s? ~ i{i + 1)e x)S* p3?dx + Qp3? { x Tx + e ) 



When the density function is rational and integral, the above method of 
solution by infinite series applies with slight modification to the case of in- 
compressibility, the equations expressing the solution being in fact much 
simplified. 

Assume 

(66) y' = i,K n x m , 

m—i 

and equate to every term in the equations which contains C m unless multiplied 
by X , but 

(67) — C m = K m . 

The main equations in the general solution then take the following forms. 

3 °° / 1 hi 

2i + l£^i K ' \m — i + 2 m — t + 3 

(68) 

h \ , gt - 2 j_ 

+ m _i + 4 +--- / /+2i + l^ i - 2 2k" 

Lim = 3 ( a A m — 2 + 4K1 Am— 3 + ~%K<l -ttm— 4 H~ ' * " ) 

(69) 3 

— -, tt-t ; — r— -T ( h\ Am-3 + 2k 2 Am-l + 3& 3 -4m_5 + * * ' ) , 

{m — i) (m + i -\- \ ) 
mK m + i(i + l)D m = bk' 2 \ mLm + him — l)Lm-i + ki(m — 2)L m - i 

+ • • • — 77 ( ^4m_2 + &1 Am-3 + h ^m-4 + ■ • • ) , 

(70) r k J 
K m + (m + 1 )D m = 6&' 2 I X m + fti i^i + & 2 iw_ 2 + • • • 

n 1 

— p ( 5m-2 + &1 Bm-3 + &2 Bm-i + '") I > 

(71) JT< + (» + 1 ) 5,- = fcF 2 (i t - ^«) , 
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(m — i) (m + i + l)D m = bk' 2 fa Lm-i + 2fa im-2 + 3& 3 L m - 3 

+ ■■■ -f p ( Am-2 — mBm^.2 + fa (^m-3 ~ mBm-3 ) + •••) , 

(72) 2 r 

K m = — ( m + 1 ) D m + bk' 2 \ Lm + fa Lm-i + fa Lm-z + fa Zm-3 



+ * * * — T} ( ^m-2 + ^1 Bm-» + * - ' ) 



These correspond respectively to equations (55), (56), (58), (59), and (60). 

Proceeding as before, equations (72) serve to determine every K m and D„ 
in terms of Ki, Di, and Ai-i. The three equations for determining these 
(corresponding to (61)) are 



K i+ (i + l)D i = bk-[-^ 1 t( m -i + 2)(^ r2 + 
m - i + 3 + m-i + 4: + ' ' ' ) Am + \2l+l ~ ik') A ^ 2 ~ 2k' \' 



(73) h 



T,K m + 2 E (m + l)A m = 0, n E (mB m + A m ) = 0. 

m=ii m=i—% m—i— 2 

3. Case of fluidity or zero rigidity. If ju == 0, both X/m and b become in- 
finite, but their ratio is finite. Hence we multiply equations (57) by /x/X 
and introduce the notation 

(74) *_£_«*?. 

The first members of these equations become 

(75) I>C m z"' and £C m a; m , 

while the second members are unchanged except by the substitution of b' 
for b. Equations (58) therefore become 

mC m = V k' 2 mLm + fa(m — 1 )im_i + fa (m — 2)Zm_ 2 + • • • 

— ( Cm-2 + fa Cm-3 + fa Cto_4 + • • • ) 

( ?6 ) ~ J, Mm~2 + faA^ + faA m -4+ ■■•)], 

C m = b' k' I i m + &i im_i + &2 im-2 -f- • • • 

TC 1 

— p ( -Bm-2 + &1 5m_3 + &2 jBm-4 + ) J , 
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which must be satisfied for m = i , i + 1 , • • • . As before, these equations 
are identical when m = i, giving 

(77) C.-V*{u-*&) 

corresponding to (59) ; while for other values of m they must be solved simul- 
taneously. The simplest procedure seems to be to eliminate C m as before, 
giving an equation similar to the first of (60) but with the first member zero. 
If the solution is carried out with successive values m = i + I, i + 2, •••, 
having due regard for equations (50), (51), and (52), every coefficient may be 
determined in terms of those of lower order, and ultimately in terms of those 
of lowest order. It is found, however, that when (76) are satisfied for 
m = i + 1 , i + 2 , • • • , only two coefficients remain arbitrary; all coefficients 
may in fact be determined in terms of Z>,- and A±.i so as to reduce to zero 
every term in the equations corresponding to (57) except those of degree i . 
In order to cause these to vanish, and also to satisfy the boundary conditions, 
it is necessary and sufficient to so determine Z>,- and -4»_2 as to satisfy the 
equations corresponding to (61). Since n = the third of these is satisfied 
independently of Z)»- and Ai-z , while the first and second take the forms 



d = ft'^f^-j E (C, ~(m-i + 2)A m ) 



(78) 



x f i + h + *! + 



) 

, 2i-2 £ 1, 1 

+ 2t + 1 Ai ~* ~ 2¥ ~ ik' Ai ~ 2 J ' 

Ec m = o. 



Special case of fluid equilibrium. If n = and ii + 0, the above solution 
fails unless terms of degree i — 1 are introduced in the four series (47). This 
case is best treated independently, especially as it is possible thereby to show 
that the above general theory includes as a special case the ordinary theory of 
the equilibrium of a fluid sphere of variable density. 

It is easy to show that, if n = , equations (44) cannot be satisfied with an 
arbitrary law of density except by making 

(79) y = 0, R = 0. 

Substituting the value of R given by (45) with y = in the equation 



«> a(*sr) -«•• + "*-•. 
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and noting equation (65), there results 

i(i + l)e — x 2 -J-; , 

(81) — dX 2pX 



de 

e+X dx- 



I px 2 dx 
Jo 



a well-known equation in the theory of equilibrium of a fluid sphere of variable 
density. 

The solution by infinite series when the density function is rational and 
integral is easily carried out if terms of degree i — 1 are introduced in the 
series.* The analysis in fact simplifies greatly from the fact that 

(82) C m = and L m = 

for every value of m . 

IV. Solution of equations for cases of variable elastic moduli 

1. The differential equations and surface conditions. The differential 
equations (37) and the surface conditions (39) cover the general case in which 
X and n are any functions of r . When X and fi are constants (37) reduce to 
(44). Restoring the terms involving derivatives of X and n, there must be 
added to the first members of (44) the terms 

,„„ s d\ „ d(ex)du ( da\du 

^ *«dx + 2 *^r£> x { e + x d^)fx- 

The surface conditions (46) are unchanged, since they do not involve deriva- 
tives of X or n . 

2. Elastic moduli assumed to be rational integral functions of r. We now 
assume that X and M are expressible by polynomials in x: 

(84) X = X (1 -Mia: + *2a; 2 + ■■■), p = /x (l + h x + h x 2 + • • • ), 

and proceed to examine the effect of the terms (83) upon the solution by infinite 
series. The density is still assumed to be given by (40). 

The substitution of the assumed series (47) for y , z , e , a reduces the differ- 
ential equations to two equations similar to (57), but with the addition of 
terms derived from (83) and (84). If we let 

(85) b = gp m a/ no 

(which includes the definition of b given by (53)), the second members of the 
new equations will be identical with those of (57), while the first members will! 
be the following: 

* This is unnecessary if fci = in the density function. 
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EJ(^+2)mC m + i(i+l)2) m + r^m5 1 + 2(m-l)< 1 1c m _ 1 

+ r^m* 2 + 2(m - 2)felc m _,+ p%w, + 2(m - 3)< 3 JC m _3 + • • • 

+ i (t + 1 ) (d Z>^! + < 2 Z> m _ 2 + • • • ) 

+ 2[»rfi < 4 m _ 1 + 2(m- 1)< 2 J m _ 2 + 3(m-2)< 3 ^ m -3+ •••]U m , 

2 { (^ + 2) C m + (m + l)D m + (^ Sl + 2d) CL-i 

+ (^% 2 + 2< 2 ) C_« + (-*, + 2t 3 ) C„_ 3 + • • • 

+ mix D„^i + (m - 1 )< 2 D„_ 2 + (m - 2)d A»-3 + • • • 
+ *iM»-i + (m-l)B mr . l ] + 2t t [A m ^ + (m-2)B m ^,] 

+ 3< 3 U»- 3 + (m -3)£ m _ 3 ] + •••la-. 



(86) 



These expressions contain no terms of degree lower than i, except terms 
containing h as a factor. Unless t\ = the series for y , z , e , a must contain 
terms of degree i — 1 . No terms of degree less than i + 1 are thereby 
introduced into the second members of the equations corresponding to (57), 
while the terms of lowest degree in the first members are 



(87) 



\(- + 2 Vi - l)C*_i + i(i + 1)B W + 2(t - l)*i ^l**"*, 
[(^ + 2) C M + Wt-i + dU<_ 2 + (»• - 2)5^)] »*-»; 



which must be made to vanish by determining C<_i and Z>i_i in terms of ^4<_ 2 , 
remembering (50). 

In the subsequent analysis it will be assumed that t± = si = 0, — a reason- 
able assumption in any geodynamical application, since it is natural to suppose 
that d\/dx and dp/dx vanish at the center. 

Using (86) instead of the first members of (57), and assuming the new equa- 
tions to be satisfied identically, we obtain a pair of equations whose second 
members are identical with those of (58) and whose first members are the 
following: 



(88) 
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Q* + 2^ mC m + t (t + 1 )D m + \j Q ms2 + 2 (m - 2) J C m _ 2 

+ ^[«wi+2(m-3) UCU-I+ ••• 

+ i(t + l)(< 2 Z) m _ 2 + < 3 Z) m _3+ •••) 

+ 2[2(m- l)t,A m ^ + 3(m-2)t,A m ^ t + •■■}, 

(j q + 2) C m + (m + 1 )D m + (£-% 2 + 2* 2 ) C^ 2 

+ (^53 + 2^0^.3+ ••• + (m-l)* 2 A^ 2 

+ (m-2)< 3 A»_3+ ••• +2f 2 [J ro _ 2 + (m-2)5^ 2 ] 
+ 3U4^-3 + (m- 3)5^]+ ••-. 
For m = i the two equations are identical, giving 

(89) ( ^ + 2) d . + ( i + 1 ) Di + MlrJ_) * 2 ^ = &jfe* (i, - |> J,_ 2 ) . 

For m > i they may be solved for C m and Z> OT , giving these in terms of coeffi- 
cients of lower order, and thus ultimately in terms of C,-, Z),, Ai-. 2 . If the 
resulting series are convergent, the solution may be completed as in the case 
of uniform elasticity already treated; i.e., by substituting the values of C m 
and D m (expressed as linear functions of C,-, Z>,-, Ai-2) in (89) and in the 
boundary equations (46), thus obtaining three linear equations for determining 
Ci, D it Ai-2. These equations, corresponding to (61), are the following 
( Xi and m being surface values of X and fx) : 



Cl +2 ) ft 



4(i — 1 ) 
+ (i+l)D t + -H — } - h Ai-2 



+ '") + 2i+l At ~ 2 ~ 2V ~ ik' Ai - 2 \ ' 



i+2 
(90) ' 

+ h- 

^ m - 1 + 3 

-EC B + 2 E (m + l)A m = 0, mi £ (mB m + A m )=0. 

Ml i$f=i m=i— 2 m—i—2 

3. Convergency of series. When the above process of solving for C m and 
D m is carried out, it is probable that the resulting series do not converge unless 
important restrictions are imposed upon the constants s m , tm. The values 
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of C m and D m , so far as these depend upon terms occurring in (88), are as 
follows : 

( m — i)(m + i + l)D m 

= - t 2 [ 4Cm- 2 + (m 2 - m - i 2 - i)Dm-2 - 2 (m - 2) Am-2 

+ 2m(m- 2)5^2 ] 

— t 3 [ %Cm-z + (m 2 — 2m — i 2 — i)^ — 3(m — 4)^4 m _ 3 

+ 3m(m-3)5 m _ 3 ] 

— U [ SCm-t + (m 2 — 3m — i 2 — i)Dm-i — 4 (m — 6)^4 m _4 

+ 4m(m - 4)B m _ 4 ] 

(m-0(m + z + l)f^ + 2)c m 

(91) = - (m - i)(m + i + l) — (s 2 CU_ 2 + s 3 Cm- 3 + $i Cm-i + • ■ ■ ) 

Mo 

-fe[2{ (m + l)(m-2) - t(t + 1) } C^ 2 + 2»(t + 1)A-, 

+ (4(m-l)(m + l)-2»(* + l))4 1 _j 

-2»(» + l)(m- 2)fl m _ 2 ] 

-< 3 [2{(m + l)(m-3) - t(t + 1) } Ck-i + 3»(» + 1)A*-* 

+ (6(m-2)(m + l) - 3»(» + 1))4*_, 

-3i(i + l)(m-"3)5 m _ 3 J 

-< 4 [2{(m + l)(m-4) - t(t + 1) } C_i + 4»(* + 1)Z>^, 

+ (8(m- 3)(m + l) -4»(* + l))4_« 

-4t(» + l)(m-4)U^_»] 
................. . ^ 

So far as convergence is concerned, the governing terms are those of highest 
degree in m , after every A„ and B n has been replaced by its value in terms of 
C„ and D„ by equations (52). After division by (m — i) (m + i + 1 ) these 
governing terms (of degree in m) are as follows: 

D m = — <2 Dm— 2 — <3 Z),„_3 — £4 Z)„,_4 — • • • , 



(92) 



_ Xp *2 + 2jUQ <2 ^ _ Xo 33 + 2jU < 3 

\ i o l^m—2 \ i o ^m^-3 

\ + 2/tio Ao + 2ft 
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Inspection of these terms shows that it is possible to insure the convergence 
of the series C m and D m by making all the governing terms zero except one in 
each equation and making the coefficients in the non-vanishing terms less 
than unity in absolute value; for in such case each series will ultimately 
approach coincidence with a convergent geometric series. The solution is 
thus practicable if X and n are expressed by binomial functions 

X = X (l +s n x n ), 
(93) 

H = /* (1 + tnX n ), 

in which n is a positive integer and s n , t n are less than unity in absolute value. 
Even with this restriction the solution is of very considerable interest, since 
it supplies a very important generalization in the case of the geodynamical 
problem of the elastic yielding of the earth. 

V. Particular solutions 

1. Numerical applications. The general solution above given has been 
applied in the computation of a considerable number of numerical results. 
The data used in these computations have been chosen in approximate con- 
formity with the known dimensions and properties of the earth, and the results 
fall into several series designed to show separately the effects of different 
assumptions regarding density, compressibility, and rigidity. Since the labor 
increases considerably with each added term in the density formula, the com- 
putations have been restricted to the case in which the formula is a binomial. 
Fortunately this simple formula is capable of expressing a fair approximation 
to the facts for the earth so far as known. It has already been pointed out 
that the above solution for variable elasticity is restricted to the case in which 
X and n are expressed by binomials. 

While omitting numerical details, it seems desirable to show the form to 
which the main formulas have been reduced for the purpose of computation. 
These will be given in a form sufficiently inclusive to cover all the special 
cases for which computations have been made. 

2. Formulas for density and elastic moduli. The simplified formulas as- 
sumed for p , X, fi are as follows: 

(94) p-po(l-Jb*), 

(95) X = X ( 1 - hx 2 ) , n = po ( 1 - hx 2 ) , 

in which k and h are positive and less than 1. It will be observed that (95) 
fall under (93), and also that they make X//j constant. 
The surface values are 

(96) Pi = Po(l-i), 
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(97) Xi = X (l-A), m = ix (l-h), 

while the mean density is 

(98) A» = Po(l-f*), 
so that 

(99) V = SL = S 



p OT 5 — 3k ' 

The numerical values to be assigned to the constants k, h, po, X , Mo in 
order to represent as nearly as possible the facts for the earth will be considered 
later. 

3. Disturbing potential of second degree. The particular solutions which 
will be considered all fall under the case i = 2 , which covers the actual case 
of tidal and centrifugal forces. The lunar or solar tidal potential is in fact 
expressed by (29) with 

(100) S t = S 2 = cos 2 0' - J, 

in which 6' is the zenith distance of the disturbing body. Referring to (42) 
and (21) it is seen that, with this value of S2, e denotes the ellipticity of a 
surface which in the unstrained body would be spherical and of radius r, 
while a denotes the angular displacement of a radius vector drawn to a particle 
for which $' = 45° . 
In the static problem we may take 6' = 8. 

4. Differential equations and boundary conditions. We now return to 
equations (37), (38), and (39), and note the form taken by their solution in the 
case represented by (96), (97), and (100). After the introduction of x as 
independent variable the differential equations are 

(ioi) = ^( x ^-?X A***" «*»«*). 

/•v , o n 1 djxz) ( da\dn gpa 

(\ + 2 fl )y + n-^ + x{e + x-^)^= — (R-np m < X x>), 

(1 

*=![kH-<-^)-°f^-iH 

xj **** 2"i 

(103) \y + 2 li ^~^ = 0, m(* + *|:) = 0, (whenx = l). 



in which 

(102) 

+ 
while the boundary conditions are 
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These are obtained from (44), (45), and (46) by putting i = 2 and adding 
the terms (83) to the first members of (44). 
Equations (23), defining y and z, now become 

1 d(ex?) Idiots?) 

5. Solution of equations. In the present case only even powers of x are 
required in the assumed series (47), which become 

y = C 2 3? + C 4 x 4 + • • • z = D 2 s? + Di x 4 + • • ■ 

(105) 

e = A + A 2 x 2 + A t x 4 + • • • a = -Bo + B 2 x 2 + B t x 4 + • • • . 

Because of (104) the coefficients in (105) are related as follows: 

(106) C m = (m + 3)A m -6B m , D m = - A m + (m + 2)B m . 

This requires that 

(107) B = iA , 

while f or m = 2 , 4 , • • • , A m and B m may be expressed in terms of C m and D m 
by the formulas 

,,„„, . (m + 2)C m + 6D m C m + (m, + 3)D m 

(1U») A m = —TZ, 7" M » -Dm = 



m(m + 5) ' m m(m + 5) 

which are equations (52) with i = 2 . 

The solution now proceeds as outlined in equations (54)-(61), but taking 
account of the terms which must be added when X and n are variable. These 
terms are obtained from the results given in Part IV, noting that all the con- 
stants s m , tm in (84) vanish except s 2 , t 2 , each of which is now replaced by — h . 

The value of R becomes 



(109) R = p (i2X 2 + i 4 x 4 + •••), 

in which 

m m + 2/~ 1 ~5^° 2k" 



(110) L 2 = | E (0» - mA.) (^ - ~rr, 1 h i^ - 



while f or m = 4 , 6 , • • • 

3 3 

(111) h m = y— . . : ;rr ( O m — 2 «Cm_4 Jrc^ro-^) +^€m— 2 r^nM- 

(m — z) (m + 6) 5 

These correspond to equations (54), (55), and '(56). 

The substitution of (105) in (101) gives the equations corresponding to 
(57); the second members are in fact identical with those of (57), while the 
first members are given by (86), proper substitutions being made for h\, k 2 , 
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• • • ,s\,Hi, • ■ • , h , h , •'• , and b being defined by (85). These equations are 
£ [ ( h + 2 ) mC m + QD m - ^ w^+2(m- 2) ) _*&_, 

- 6^2)^-2 - 4 (m - 1 ) ^^2 z m 
= 6&' 2 ( 1 - kx 2 ) £ |"mi m a; m - 3<7 m Q x m + 2 - ^ib; m+4 ) 

(112) -£4. *■*•], 
Z [Qjj + 2) 0, + (to + 1 )D m - (£-° + 2) AC m _ 2 

- (m- l)AZ) m _2-2/j(^ m _2 + (m-2)5 m _ 2 )Ja;' B 

= bk' 2 (1 - fa 2 ) Z [i m * m - p&. * m+2 ] . 

the summations extending to all values of m which are consistent with (105). 
Both equations thus take the form of series of even powers of x, beginning 
with x 2 . To satisfy them identically it is necessary to satisfy the following 
equations for m = 2 , 4 , 6 , • • • : 

+ 6Z) m _ 2 + 4(m- 1)4-«J 

r 8 3 

= bk' 2 mL m — k ( to — 2 ) Z m _ 2 — C m -z + -z W m -i — ■= k 2 C m -t 

(113) - jp (Ar^-2 - kAm^i) \, 

Q* + 2) C m + (ffi + 1 )D m - h [(^ + 2) 0^2 + (m - 1 )B—« 

+ 2^ m _ 2 + 2(m- 2)£ m _ 2 J 

= &&' I i m — kLm-2 — T> ( 5m_2 — kBm^i ) . 

For to = 2 the two equations are identical, giving 
(114) Q? + 2) C, + 3Z> 2 - 2hA = 6F 2 (z 2 - ^>4o) ; 

while for to = 4, 6, • • • , they may be solved for C m and Z> m , giving 
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bk' 2 T 8 3 

(m — 2) (to + 3) L 5 5 

+ p ( ^4 m _ 2 - mjBm-2 — & ( Am-i — raB m _ 4 ) ) 

,/ n 4C m _ 2 - 4Z) CT _2 \ 

(115) +A V 2 (^-2)(m + 3)/' 

(m"o + S ) Cm = ~ (m + 1 )Dm + 6fc ' 2 \ Lm ~ kLm ~ 2 

- p (5 m _ 2 - **„_«)] + A [Qj + i) C_ 

4(m - 1)0^2 + (m 3 + 2m 2 -9TO + 14)Z> ro - 2 "l 
+ (m-2)(m + 3) J* 

It should be noted that, in the terms containing h , Am- 2 , and Bm- 2 have been 
replaced by Cm- 2 and Z>m_2 by means of (108). In applying these equations, 
L m and Lm-i are to be replaced by their values as given by (111), except in 
the case of L 2 which is to be eliminated by means of (114). In this way every 
C m and D m may be computed in terms of coefficients of lower order, and 
ultimately in terms of C 2 , D 2 , Ao . 

For purposes of computation the formulas for C m and D m may conveniently 
be expressed explicitly in terms of like coefficients of the three preceding orders. 
This is possible when m > 6; when m = 4 and 6 special formulas are required. 
Convenient working formulas may be written as follows: 

*-£(ft-JftH[(£+»)a + <A-2u.] 

-y[(^ + 2)<7 2 + 3Z> 2 -2L4 ] + a[(^ + 2) C 2 - ^(C 2 - 3Z) 2 )] . 

hJr' 2 r If 1 2k 2 

Dt = ~k L Ci " S ( 61C2 + 30Z>2 } + it Ao 

- jjp(z) 4 -*(C, + 12D 2 ))] +J(C 4 + &D 4 ), 

(H7) . . /,jfe' 2 r JS- 

^ + 2 J C 6 = ^- [ - 4C 4 + 6Z) 4 + 35 (286C 2 - 564fl 2 ) 
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Form = 8, 10, ••• , 

D m = bk' 2 \-. ^ — x - k; 

L(m - 



2(4m 2 - 7m - 41 ) CW_ 4 + 602)^4 
2)(m + 3) /c 5(m-4)(m-2)(m + l)(m + 3) 



3 (m 4 - 8m 3 - m 2 + 76m + 28) C— «+ 36 (m 2 - 3m - 4)Z) ro -« 
+ 5(m-6)(m-4)(m-2)(m-l)(m + l)(m + 3) 



k' 



% ( D„^2 20^-4 + (m 2 -m -6)A^ 4 \1 

fc'\(m-2)(m + 3) (m -4) (m - 2) (m+1) (m+3)7 J 

/ 4(7^-42)^ \ 

+ ft V m - 2 + (m-2)(m + 3)y' 

(118) ^°+2)c ra = M'{p%^; 

v \Mo / L(™ -2)(m + 3) 



4(8m 2 - 8m - 97)C m - i - 12(4m 2 - m - 44)Z) m _ 4 
+ 5(m-4)(m-2)(m + l)(m + 3) 



-k 2 



12 ( m 4 - 6m 3 - 14m 2 + 69m + 46 ) CU-e 

- 18 (m 2 - m - 8) (m 2 - 3m - 14) Z> OT _ 6 



5(m-6)(m-4)(m-2)(m-l)(m + l)(m + 3) 

n( C^ 2 , (m 2 -m -8)0^ + 121)^-4 \"| 

&'V(m-2)(m + 3) (m -4) (m - 2) (m + 1) (m+3),/J 

It is seen that, in (118), the governing terms as regards convergency are 
those of degree in m. Thus the formulas ultimately approach the forms 

ai9) D m = fcIV-2, C m = hC^i, 

which indicate convergency if h < 1 ; this of course includes the case of uniform 
elasticity (h = ) . 

It remains to satisfy (114) and the boundary conditions (103), L^ being 
replaced by its value (110). These reduce to the following three equations, 
which, after substitution of the values of the summations, become linear 
equations in d, Z) 2 , Aq\ 

(- + 2) C 2 + 3Z) 2 - 2hA 

~ 6 * L25V4^ m "t B + 5 j 
/torn 3,/^ q»-2Z) CT ^ C m -2D m \ , /2- n\ c 1 
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/Xo x .\f „ , 4^ C m + 3D m 24 - C m -2D m 
V Mo / a 5 2 hi 5 2 to + 5 

2 O 2 TO O2TO+O 

The cases for which numerical results have been computed are all covered 
by the foregoing equations; the specific formulas for each case being easily 
written by assigning particular values to k , h , and n . 

6. Case of uniform density and elasticity. For the purpose of comparison 
with previously known solutions, it is of interest to show the form assumed 
by the main working formulas in the case h = , k = , both for unrestricted 
compressibility and f or X = °o . 

(A) Compressibility unrestricted. Putting h = 0,k = 0,k' = l, equations 
(118) become 

Dm = (to-2)(to + 3) (<7m - 2 " nZ)m - 2) ' 
(121) 

Q + 2 ) C '" = (m-2) 6 (TO + 3) (-( 4 + ^ C - + 6Z) - ) ' 

and hold f or to = 4 , 6 , 8 , • • • . Successive applications of these equations 
determine every C m and D m in terms of C2 and D 2 . These values are to be 
substituted in (120), of which the second and third are unchanged while the 
first simplifies as follows: 



(122) 



\/x / |_25\lr to 2 to + 5 / 



+ 



(f-l)^-l] 



(B) Incompressibility. When X = » we proceed as indicated in the general 
solution, equations (63)-(69). Equations (121) then become 

(123) ft.-- ,.. *ffi?, QW K m = 6W ^ 2 



(to-2)(to + 3)' m (to-2)(to + 3)' 

while (120) reduce to the following: 
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Static problem. When n = equations (123) and (124) result in the well- 
known solution for the static strain of a homogeneous incompressible sphere. 
Thus equations (123) show that, if n = 0, every D m vanishes except D 2 , 
that every K m vanishes except K 2 and K* , and that 



2 • 



(125) Ki = f bD 
Equations (124) are thus reduced to the following: 

K 2 + (3 + ^&) D 2 - jbAo + \bc = 0, 

(126) A- 2 + (y + f&)#2 + 2^o = 0, 

from which 

^o_ 86 D 2 = lb K 2 _{2 + Zb)b 

{ ' c 38 + 46* c 38 + 46' c 38 + 46^ 

Also, from (107) and (108), 

A 2 3 D 2 36 



(128) 



c 7 c 38 + 46' 

B 2 5 D 2 56 



c 14 c 2(38 + 46)' 

Bo A 46 



c 2c 38 + 46 ' 
so that the final solution is 

e 8 - Sx 2 a 8 - 5x> 



(129) 



c « c 

38— + 4 
gpa 



^38 — +4) 



the known solution for the static problem. 

Oscillation of an incompressible fluid sphere. It is of interest to note the 
simple form assumed by the solution in the case p = constant, X = » , 
H = 0, n =j= 0. This is best treated by direct use of equations (44), (45), and 
(46), noting that y = but \y 4= 0, and that the second of (46) is satisfied 
by the condition p = 0. The simplified equations are the following: 

d(\y) ( dR 



-(.£-*-). 



(130) da; 

Xy = ga(R — wpcKr 2 ) , 
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(131) 




R = p{e- 


fei - fc)* 2 , 


(132) 




\y = when x = 1 


Equations 
(133) 


(130) give 


diax 2 ) 
dx 


- ex 2 = 0. 


But since y 


= the first of (48) gives 




n!U\ 






d(ex*) 



v *^, ""- - dx • 

Rejecting negative powers of x as irrelevant to the physical problem, equations 
(133) and (134) give 

(135) e = 2a = constant = e\. 
Hence 

(136) R = p(ie 1 -$c)x 2 , 

(137) Xy = ffpa((|-|)ei-|cjar 2 . 
Finally, from the surface condition (132), 

(138) •»&..*. 
v c c 4 — 5ra 

Numerical applications of this equation will be given later. 

7. Quantities measuring the strain. Two quantities are of especial im- 
portance because their values enter into the observational data from which 
the actual strain of the earth is inferred. One of these is the surface ellipticity 
e\, while the other depends upon the change in the principal moments of 
inertia caused by the strain. 

Surface ellipticity. This may be computed by the formula 

„„ m v 1 a a i 2 y> "m + oD m 3 ^-> C m — 2D m 

(139) ei = 2- A m = A Q + r 2- ~ 1- ¥ 2^ - , g— , 

obtained by the use of (108). 

Change of principal moments of inertia. Let C, A, A be the principal 
moments of inertia of the strained body (two of them being equal since Sj 
has an axis of symmetry). Then it may be shown that 

„,,« , ^ 8tto 6 r 1 ( . d(ex*)\ J 

(140) A-C^-^l P (yx*--^)dx, 

which, by the use of (104), may be written 

(141) A - C = ^- 6 f p(e + 3a)cc 4 <£r. 
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Let I denote the diametral moment of inertia of the unstrained body, and let 

(142) / = ^. 

Then 



(143) tJ p ~3~ J P^ dx > 

2 I p(e + 3a)x 4 dx 
Jo 



so that 



(144) / = 



•I 



per 4 dx 



Using the value of p given by (94), and substituting the values of e and a from 
(105), the integrals in (144) become 

(145) £ P(e + 2«)x*dx = P0 ±y m + 3B m )(^ r5 - J^), 

(146) J px*dx = p ( 5-7V 

Convenience in the computations is gained by replacing A m and B m by C m 
and D m by means of (107) and (108), which give 

(147) Ao + 3B =^A , A m + 3B m = m Z - ( m = 2,4,...), 

& lib 

so that (144) reduces to the form 

t _ a . ? V ^ + 3Z? » 14 A C m + 3Z) OT 



52 m 5(7 — 5k) 2 m + 5 

(148) 

j 2 ^ C m -f- 3X/ m 



7 — 5k 2 m + 7 

8. Deflection of apparent gravity. A quantity of importance in the applica- 
tion to the earth is the deflection of apparent gravity relative to the earth's 
surface due to the action of the disturbing forces. 

The equation of the boundary surface of the strained body is 

(149) r = a(l+«iS,), 
while that of a level surface may be written 

(150) r = a ( 1 + e S 2 ) . 
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The angle between the normals to the two surfaces (149) and (150) is 
(151) (e-e,)-^. 

The value of e may be found from the condition that the sum of the gravita- 
tion potential of the strained body and the disturbing potential is constant 
at the surface (150). 

The change in the external gravitation potential due to the strain may be 
computed by the same method as that employed in finding the internal 
potential U; it is in fact expressed by the first integral in (35) if the upper 
limit r is replaced by a . Comparing with (140) it is seen that the change in the 
external potential is equivalent to 

n *» Sy(A-C)S 2 _ Sga> IfS 2 

llMJ 2r» ~ 2r 3 M • 

The external potential of the unstrained body is 

(153) yM = 9jt 

r r 

and the potential of the disturbing forces is 

(154) W = f a , * S2 - 

The sum of these three parts, at the surface (150), to the first order of small 
quantities, is 

(155) , a [ 1+ (^- + £_ e )s 2 ], 

and this is constant if 

nam , ° -l 3// 

(156) e= 2 + 2Ma>- 

This reduces (151) to 

which is the angular deflection of apparent gravity with respect to the- 
boundary surface. 

Let q denote the ratio of this quantity to the value it would have if the- 
body were absolutely rigid; then since the latter value is %cdSt/dd , 
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VI. Application to the earth 

1. Values of constants in case of the earth. In the case of the earth the 
mean radius and mean density are quite accurately known. Expressed in 
C. G. S. units the values here used are* 



a = 6.371 X 10 8 , 
From these result the values 



5.527, y = 6.6576 X 1(T 8 . 



g == 982 , g Pm a = 3.458 X 10 12 . 

As regards the variation of density within the earth, it is known that the 
•surface density is about half the mean density, and that the diametral moment 
of inertia is about | of the value for a homogeneous earth of density p m ; i.e., 



(159) 



2ra 



£ 



Pi = 2Pm > 



I = iMa 2 . 



uft\ 


\ 

V 




>\ 



.6 
X 

Fig. 1. 



1.0 



These equations can be satisfied only approximately 
by a density function of two terms, but the ap- 
proximation can be sufficiently close to serve the 
present purpose. The formula which has been used 
is expressed by (94) with k = § . This gives a close 
approximation to the correct value of the moment 
of inertia, but makes the surface density less than 
that of actual surface rocks. A comparison with 
the well-known formula known as Laplace's law of 
density is givenf in Table I and in Figure 1. 



TABLE L Comparison of Assumed Law of Density with Laplace's Law 

U) P = po(1 - hx*) = 2 Pm (l - fx 2 ), 
D sinicx .. sin 2.46 1s 



r 


pI 


Pm 






a 


(A) 


(B) 





2.00 


1.95 


.2 


1.93 


1.87 


.4 


1.73 


1.65 


.6 


1.40 


1.31 


.8 


.93 


.91 


1.0 


.33 


.50 



* The values of p„ and y are those found by C. V. Boys, Nature, vol. 50, p. 419. 

t The formula of Laplace satisfies both of equations (159) very closely. The same is true 

of the trinomial formula 

/. 14 , , 21 .\ 

p = '"\ 1 --m x '+65* )' 

and numerical solutions using this formula are entirely practicable, but considerable additional 
labor would be imposed by the addition of a third term to the density formula. 
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The values of /* and X have been determined experimentally for many speci- 
mens of surface rocks. Concerning the values throughout the interior of 
the earth the only evidence is the comparison of the actual yielding of the 
earth to disturbing forces with the computed yielding. Computations based 
on the assumption of uniform elasticity throughout the body have made it 
quite certain that the material of the earth is on the average far more rigid 
than the surface rocks. The computations which have now been made on the 
assumption of variable elasticity are based upon formulas (95) with values 
of the constants Xo , juo , and h so taken that 

X x = m = 2.542 X 10 u C. G. S. units. 

These surface values are obtained by averaging a considerable number of 
results given by Adams and Coker.* 

2. Assumption regarding n . The value of n for the actual tidal disturbance 
is so small that the terms containing n in the formulas expressing the solution 
may be neglected in computations relating to the problem of the rigidity of 
the earth. The problem of forced and free vibrations is, however, of some 
interest in itself, and one series of numerical results will be given based upon 
a series of values of n . 

3. Summary of numerical results for static problem. The numerical results 
which have been obtained for the static problem fall into three main groups, 
the design being to show separately the influence of variable density and 
variable elasticity upon the quantities which are involved in estimates of the 
rigidity of the earth. The grouping corresponds to the following three 
assumptions : 

(1) p and n both constant; 

(2) p variable, /z constant; 

(3) p and n both variable. 

Groups (1) and (2) each consist of two series of results, one for X = °o (in- 
compressibility), the other for X = ju. The quantities of chief interest are 
Bi/c,f/c , and q , since it is these whose actual values are inferred from observa- 
tion; and the results are so arranged as to show the way in which each of 
these quantities depends upon n . 

* An investigation into the elastic constants of rocks, more especially vxith reference to cubic 
compressibility, by Frank D. Adams and Ernest G. Coker, American Journal of 
Science, vol. 172 (1906), pp. 95-123. Averages of the values found for seventeen speci- 
mens of rocks of various kinds are as follows : p = 2.542 X 10" dynes/cm 2 , X + f M = 4.218 X 10 11 
dynes/cm J . These give X//» = 0.993. 
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(1) Case of uniform density and elasticity. (Table II and Figure 2). 
Although the solution of this case is known,* the series of numerical results 
here given is more extended than has heretofore been published. Of especial 



1.60 



1.00 



.60 

















QcSi — 









.2 

H/pga 
Fig. 2. 



interest is the curve representing the results for the case X = p. It is seen 
that this curve shows an infinite discontinuity for a value of fi/gpa between 
0.12 and 0.125. This corresponds to the case of " gravitational instability " 
which is discussed by Love.f Algebraically this case results from the fact 

TABLE IL Results for the Case of Uniform Density and Elasticity 



pga 


X =m 


X = 


00 


ei 


/ 


q 


ei = / 


q 




c 


c 




c c 







1.250 


1.250 




1.250 





.05 


.732 


.815 




.847 


.322 


.10 


.247 


.554 




.641 


.487 


.12 


-2.371 


-.030 




.584 


.533 


.125 


4.955 


1.472 




.571 


.543 


.15 


.830 


.577 


.032 


.515 


.588 


.20 


.535 


.449 


.470 


.431 


.655 


.30 


.366 


.331 


.665 


.325 


.740 


.40 


.285 


.264 


.747 


.260 


.792 


.50 


.235 


.220 


.794 


.217 


.826 



* Some Problems of Geodynamics, by A. E. H. Love. See especially Chapters VII and 
VIII. It is of interest to note the following results given by Love: 



* = ^ = 1 84 
m 25 X - 84 ' 



-=- = £7 = 0.208, 
gpa 24 



5i = 0.466, 
c 



f- = 0.4275; 
c 



2041 



= 1.05, 



-£- = 0.202, 
gpa 



- = 0.522, 

c 



c 



0.436. 



X 

m 1944 

It will be seen that these conform well to the results shown in Fig. 2. In making this com- 
parison it should be ndted that Love's symbols h , k are so denned that 



c ' 



fc = 



3/ 

Mc? 



so that the quantity above designated by q is, in Love's notation, 1 + k — h . 
t Some Prqblems of Geodynamics, Chapter IX. 
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that, for certain values of the elastic constants, equations (120) may be 
satisfied with c = (i.e., with zero disturbing force) without the vanishing 
of the coefficients in the assumed developments of e and a . 



l.OO 



< .60 













r 






fX-oo) 



.1 .2 .3 .4 

Fig. 3. 

(2) Case of variable density but uniform elasticity. (Table HI and Figure S). 
The effect of the assumption of variable density is seen by comparing the 
two series of results in Table III and Figure 3 with the corresponding series 

TABLE HI. Results for the Case of Variable Density and Uniform Elasticity 
(Values in brackets obtained by extrapolation.) 



M 


X =n 


X = oo 














Pm ga 


ei 


/ 


q 


£i 


/ 


q 




c 


c 




c 


c 







.942 






.942 


.910 





.05 








.597 


.643 


.430 


.10 








.472 


.514 


.555 


.15 


.763 


.532 


.010 


.390 


.426 


.634 


.20 


.436 


.395 


.512 


.333 


.365 


.689 


.30 


.296 


.291 


.690 


.258 


.283 


.759 


.40 


.234 


.236 


.760 


[-211] 


[.231] 


[.802] 



in Table II and Figure 2. The variation of the strain with /i is similar in the 
two cases, but the strain for a given value of n is materially less in the case of 
variable than in that of uniform density. The two curves for the case X = m 
resemble each other also in the- feature of an infinite discontinuity.* 



* The numerical solution in the case of variable density involves a considerable amount 
of labor; it is, however, greatly facilitated by the use of calculating machines. In preparation 
for the solution of a series of cases with specific values of the elastic moduli, the coefficients in 
formulas (116), (117), and (118) which depend upon k and m may be computed once for all 
for values of m as high as may be deemed necessary; these are used in connection with the 
numerical values of 6 and of X//x appropriate to each specific case. Every D m and C m is thus 
computed from terms of lower order, the resulting values being linear functions of»Cj , Dj , 
and Ao . From these the summations occurring in equations (120) may be computed as linear 
functions of Cs , Di , and A*; it is then easy to complete the solution. 
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3. Case in which both density and elasticity are variable. (Table IV and 
Figure 4). A series of results for the case of variable elasticity was obtained 
by giving to the constants no and h in the formula X = /j=/*o(1 — fox?) four 
sets of values, so taken that n\ has always the value 2.542 X 10 u already cited 
as holding for surface rocks, while /*o/mi has the four values 4, 5, 6|, 10. The 
four assumed values of h are 0.75, 0.80, 0.85, 0.90. Table IV and Figure 4 
show the way in which the three quantities ei/c , f/c, and q vary with po • 



1.00 



'.60 







jO££— 








^o^ir 





Fia. 4. 

4. Rigidity of the earth. Experimental evidence as to the yielding of the 
earth to small disturbing forces involves the two quantities f/c and q . 

The value of f/c may be inferred from the earth's free nutation period as 
determined by observed variations of latitude, the prolongation of the period 

TABLE IV. Results for Variable Density and Variable Elasticity 





M0 




Mo 


ei 


/ 




h 




6 




— 




q 




Ml 




Pm ga 


c 


c 




.75 


4 


3.4025 


.294 


.605 


.440 


.217 


.80 


5 


2.7220 


.367 


.374 


.335 


.576 


.85 


6| 


2.0415 


.490 


.283 


.269 


.695 


.90 


10 


1.3610 


.735 


.211 


.207 


.779 



as compared with that for an absolutely rigid body being accounted for (as 
Newcomb first pointed out) by elastic yielding.* This line of evidence 
indicates a value of f/c between 0.27 and 0.29. 

The observed effect of lunar and solar tidal forces upon the direction of 
gravity relative to the earth at a given place should theoretically give the 

* It may be shown that the actual nutation period is to a close approximation the same 
as that of an unyielding body having the figure which the earth would assume if centrifugal 
forces were annulled. The theory leading to this result was given by the present writer in 
a paper presented to the American Mathematical Society, San Francisco Section, Dec. 20, 
1902, of which an abstract was published in the Bulletin of the Society, vol. 9 
(1903), p. 299. Similar reasoning was employed by J. Larmor in a paper in Proceedings 
of the Royal Society, London, Ser. A, vol. 82 (1909). 
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value of q. This quantity is, in fact, the ratio of the actual deflection of 
gravity to the deflection computed on the assumption that the earth is abso- 
lutely unyielding. Precise observations of this effect have been made by 
O. Hecker and others by the use of the horizontal pendulum, and by A. A. 
Michelson by measurement of the actual water tides in closed horizontal 
pipes. The observations of Hecker* gave values of q varying with azimuth, 
indicating a greater yielding of the earth in the north-south than in the east- 
west azimuth. Michelsdh's results,t however, indicate equal yielding in all 
azimuths, the value of q adopted as most probable being 0.69. 
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The values of f/c and q which are inferred from observation cannot be har- 
monized with theory on the assumption of homogeneity and incompressibility, 
and the discrepancy is only partly removed by assuming compressibility. 
The results above given show that both variable density and variable elasticity 
work in the direction of harmonizing the computed with the observed values. 
This is brought out by comparing the simultaneous values of f/c and q in 
each of the five series of results given in Tables II, III, and IV. The tabulated 
data are represented graphically in Figure 5, the five curves showing the rela- 
tion between f/c and q for each of the five series. The curves (la) and (2a) 
correspond to incompressibility, while (1), (2), and (3) are for the three cases 
in which X = fi . 

Comparing the three curves (1), (2), and (3), it is seen that if f/c lies between 

* O. Hecker, Ve rof f e nt li c hu n gen des Kgl. Preussischen Institutes, 
No. 32, Berlin (1907). 

t Michelson's preliminary observations are described in a paper entitled Preliminary 
results of ■measurements of the rigidity of the earth. Journal of Geology, vol. 22 
(1914). See also Astrophysical Journal, vol. 39 (1914). These measurements 
gave 0.71 as the best value of q. (For a correction to the results originally published see 
Science, October 3, 1919, p. 327.) The value 0.69 was obtained from a more extended series 
of measurements, the aceount of which has not at this date (October, 1919) been published. 
The writer is indebted to Professor Michelson and Professor Henry G. Gale for the com- 
munication of the general result in advance of publication. 



42 L. M. HOSKINS [January 

0.27 and 0.29 (as indicated by latitude observations), q should fall within the 
limits 

0.74 and 0.72 in case (1), 0.72 and 0.69 in case (2), 0.69 and 0.66 in case (3). 

Materially greater values of q are indicated by the curves (la) and (2a). 

Considering all the computations that have been made, the best agreement 
with the facts of observation results from the assumptions 

p = Po ( 1 - f x 2 ) , X = p = p ( 1 - .85a; 2 ) , 

• 

the value of po being based upon the known value of p m , and that of p being 
such as to make pi agree with the value found experimentally for surface rocks. 
This surface value is 

Px = 2.542 X 10 11 dynes/cm 2 , 

giving for the value at the center 

po = pi/.15 = 1.695 X 10 12 dynes/cm 2 . 

5. Forced oscillation. The numerical results given above are for the static 
problem, all terms in the formulas which contain n being regarded as negligible. 
Inspection of equations (118) shows, however, that the same routine of com- 
putation can be followed when n has any known value. It is thus possible 
to determine the forced oscillation due to a disturbing potential of any known 
period. 

6. Free oscillation. If an oscillation of the type assumed in the solution 
is possible in the absence of disturbing forces, equations (120) must be satisfied 
with c = . Obviously this is possible only for particular values of n , and 
the only method which suggests itself of determining such particular values 
is by interpolation in a series of results computed with assumed values of n . 
This method, though somewhat laborious, is entirely feasible, and has actually 
been carried out in one case with results which will be given. 

If we consider a continuous series of forced oscillations with the same 
intensity of disturbing force but with period decreasing from a> (corresponding 
to static strain), the period of possible free oscillation will be that at which 
the forced oscillation reverses its phase while its amplitude passes through » . 
If, therefore, for each of a series of values of n , beginning with , the values 
of d/c , Dz/c , Aq/c are computed, each of these will increase in magnitude 
with increasing n, approaching positive or negative infinity as n approaches 
the value corresponding to the free oscillation, then changing sign and de- 
creasing in magnitude. The required critical value of n may be found by 
interpolation in one of these series, preferably by taking reciprocals and 
interpolating for the zero- value. 
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The case for which computations have been made is that in which the 
density and elastic moduli are represented by equations (94) and (95) with 
k = f and h = 0.85. Results for n = have already been given in Table 
IV, and further computations have been made for n = 1 , 2 , and 3 , giving 
the results shown in Table V. 

TABLE V. Results for Case of Forced Oscillation 
[Results in brackets obtained by interpolation for case of free vibration.] 



n 


Period in 
seconds. 


C 2 
c 


Z>2 

c 


c 


ei 
c 


c 

«i 



1 
2 
[ 2.77 ] 
3 


5060 

3578 

[3040] 

2922 


.795 
1.668 
6.820 

[± • 1 
-44.04 


- .540 .381 

- 1.146 .608 

- 4.687 1.492 

[ ± 00 ][ ± oo ] 

+ 30.90 - 5.857 


.283 

.436 

.948 

[± » ] 

-2.93 


3.53 
2.29 
1.06 
[0] 
-.34 



The interpolation for the critical value of n may be accomplished to a good 
approximation by means of the graph (A) shown in Figure 6, which represents 
c/ei as a function of n; the zero-value falling at about n = 2.77 . 

From the meaning of p is given by (27) , it is obvious that the period of a 
complete oscillation is 27r/p or 2ir ^la/ng . Putting a = 6.371 X 10 8 cm., 
g = 982 cm/sec 2 , the period in seconds is 5060/ Vn. The period corresponding 
tow = 2.77 is thus about 3040 seconds. 
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It is of interest to compare this result with that given by equation (138) 
for a homogeneous incompressible fluid sphere. The critical value of n for 
this case is 0.8, corresponding to a free oscillation period of 5650 seconds. 
The graph of c/ei as a function of n is the straight line (B) shown in Figure 6. 



